Two Kinds of the Coexistent States in One-Dimensional Quarter -Filled Systems 

under Magnetic Fields 
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The coexistent state of the spin density wave (SDW) and the charge density wave (CDW) in the 
one-dimensional quarter-filled system and with the Coulomb interaction up to the next-nearest sites 
under magnetic fields is studied, ft is found that, in the coexistent state of 2fcp-SDW and 2/cp-CDW, 
the charge order is suppressed and 2fcp-SDW changes to 2fc F -SDW having the different alignment 
of spin under high magnetic fields, whereas, in the coexistent state of 2fep-SDW and 4te-CDW, 
4/cf-CDW still remains and 2fcp-SDW is suppressed. The critical temperature of the charge order 
is higher than that of the spin order when the inter-site Coulomb interaction is strong. These will 
be observed in experiments such as the X-ray scattering measurement and NMR. 

PACS numbers: 71.45.Lr, 75.30. Fv, 75.50.Ee 



I. INTRODUCTION 

Organic conductors such as (TMTSF) 2 X and 
(TMTTF) 2 X (X=C10 4) PF 6 , AsF 6 , Re0 4 , Br, SCN, 
etc.) are the quasi-one dimensional quarter-filled systems 
and exhibit various kinds of ground states, for example, 
spin-Peierls, spin density wave (SDW) and supercon- 
ductivity. dj| In (TMTSF) 2 PF 6 , the incommensurate 
SDW occurs at T = 12 K, where the wave vector is de- 
termined by NMR experiments f||| as (0.5, 0.24, 0.06). 
Pouget and Ravy observed the coexistence of 2fc F -SDW 
and 2fcp-charge density wave (CDW) by the X-ray scat- 
tering measurement, || where A; F = 7r/4a is Fermi wave 
number and a is the lattice constant. The critical tem- 
perature (Tcdw) a * wmcrl 2/cp-CDW is observed is the 
same temperature as Tsdw- 

On the other hand, it is found that SDW in 
(TMTTF) 2 A" (X=Br and SCN) has the commensurate 
wave vector, (0.5, 0.25,0), as observed in the measure- 
ments of 13 C-NMR @ and 1 H-NMR §. From the angle 
depenence of satelite peak positions of 1 H-NMR |@jg, 
the alignment of the spin moment along the conductive 
axis (a-axis) is obtained to be (|, 0, L 0), where the arrow 
means the spin moment. In (TMTTF) 2 Br, 4fc F -CDW 
accompanied by 2fc F -SDW is found in X-ray scattering 
measurments, where Tcdw ~ 100 K and Tsdw ~ 13 K. 
1 

In the one-dimensional quarter-filled systems, the 
CDW-SDW coexistent state is shown to be caused by 
the interplay between the on-site Coulomb interaction 
(U) and the inter-site Coulomb interaction (V). |9j-|l3j| 
The inter-site Coulomb interaction plays important role 
for the charge order. Mila |L4| has estimated U/t ~ 5 
and V/t ~ 2, where t is a transfer integral. Seo and 
Fukuyama || showed that the ground state becomes the 
coexistent state of 2fc F -SDW and 4fc F -CDW in the one- 
dimensional extended Hubbard model with V. It is found 
by Kobayashi et al. (l(]|llj an d Mazumdar et al. Jl^] 
that the coexistent state of 2fc F -SDW and 2fc F -CDW is 
stabilized when the next nearest neighbor Coulomb in- 



teraction (V^) and the dimerization of the energy band 
are considered. Thus, two kinds of the coexistent states 
(2fc F -SDW-2fc F -CDW and 2fc F -SDW-4fc F -CDW), which 
are observed by the X-ray scattering mesauerments in 
(TMTSF) 2 PF 6 and (TMTTF) 2 Br, can be explained by 
Seo and Fukuyama g] and Kobayashi et al. |10|,[ll| and 
Mazumdar et al. respectively. 

Recently, we have studied the effects of the magnetic 
field (H) on the coexistent state of 2fc F -SDW and 4fc F - 
CDW for both cases of strong and weak coupling of the 
correlation between electrons, where the one-dimensional 
quarter-filled extended Hubbard model with V is used. 
p8| It is found that although the spin order is sup- 
pressed at high fields, the charge order still remains in 
the strongly coupling systems (U/t ~ 5). When the cou- 
pling is small (U/t ~ 1.5), both orderings of 2fc F -SDW 
and 4/c F -CDW disappear at the critical magnetic field, 
which is the same as the Pauli par amagnetic field in the 
spin-singlet superconductivity. |1^ [TtJ 

In this paper, we examine how the coexistent state 
of 2/cp-SDW and 2fc F -CDW and the one of 2/t F -SDW 
and 4/c F -CDW are affected by magnetic fields and tem- 
peratures. We use the one-dimensional extened Hub- 
bard model with V and V^, where the dimerization is 
neglected for simplicity, because the coexistent states 
are stabilized without the dimerization. We use param- 
eters, U/t = 5.0, which is indicated in (TMTTF) 2 X 
((TMTSF) 2 X). (HUD In most cases we set V/U < 1 
and V 2 < V. We study the case when the anisotorpy of 
the spin is strong, because it is found that the easy axis of 
the spin in the quasi-one-dimensional organic conductors 
is 6-axis. |2f| In order to clarify the ground state under 
fields at finite temperatures, we show the V-H, V2-H, 
V-T and V2-T phase diagrams, which enables to discuss 
the effect of pressure and the ordering temperatures of 
the SDW and CDW. It is expected that the Hubbard 
model with large U/t ~ 5 is qualitatively understood by 
the Ising model with the perpendicular field, as we will 
show below. 

The _ff-dependence of the antiferromagnetic state in 
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one-dimensional Ising model has been studied in the 
mean field approximation, where the componet of the 
spin along the c-axis is considered. |2^j2^] When the mag- 
netic fields is applied perpendicular to the spin moment, 



S z (H,j)/S z (0,j) = VHW- 



where S z (H,j) is the amplitude of the spin moment at j 
site at H — and H® is the critical field at which the or- 
dering of the antifcrromagnetic state disappear (H® oc J, 
where J is the coupling constant). 



II. FORMULATION 

We study the one-dimensional extended Hubbard 
model, 



H = JC+U + V + V 2 , 

t = -t^(4 )<r Cf + i )(T + h.c.) 



i, a 



V 2 = F 2 E 



(2) 

(3) 
(4) 
(5) 
(6) 



where /ie = eU/2moc is the Bohr magneton and cj a is the 
creation operator of a spin electron at i site, ni, ff is the 
number operator, g = 2, i = 1, • • • , N$, iVg is the number 
of the total sites and a =j and j. The notation in this 
paper follows Seo and Fukuyama. || The magnetic field 
is applied to the x-axis and {H x ) aa i — H(a x ) aa >, where 
H is the strength of the magnetic field and a x is Pauli 
spin matrix. We consider the quarter-filled case. 

The interaction terms, IA, V and V2 are treated in the 
mean field approximation as 

^ M = EEi' 5 TT(Q^)^ t (^ I l)c(fc :c -g,i) 

+ pl l (Q,T)tf(k x -Q,T)C(k x ,m 
V M = (jj) E y Ee'* Qa {P™(Q>T)&(k x ,<j')C(k c -Q ) 

+ p*, a ,(Q,T)C\k x ,a)C(k x -Q,a)} 

k x ,(T,a' Q 

+ P :, a ,(Q,T)C\k x ,a)C(k x - Q,a)} 

where I = U/N$. The self-consistent equation for the 
order parameter p aa (Q,T) at finite temperature, T, is 
given by 



r (Q, T) = lJ2 < C\k Xl a)C(k x - Q, a) > . (10) 



In eqs. (7)-(10), we take the order parameters, p aa (Q, T), 
only between electrons with same spins. We do not con- 
(1) sider the case of the mean field, p aS (Q,T) = IJ2k < 



C^(k x ,a)C(k x — Q,a) > with a ^ a. We neglect the 
effect of the spin tilting to the x — y plane by setting 
Pas{Q,T) = in this mean field approximation. This 
simplification corresponds to the assumption that the ro- 
tational symmetry in the spin-space is broken and that 
the z-axis is the easy axis although the spin in this Hub- 
bard model is isotropic. The magnetic field is applied 
perpendicular to the easy axis when H =/= 0. 

In the Fulde-Ferrel-Larkin-Ovchinnikov state in the su- 
perconducutivity p6| , ^7f , the total moment of the Cooper 
pair is changed by the magnetic field. Similar situa- 
tion may occur in SDW and CDW, i.e., the wave vec- 
tor Q may be changed by the magnetic field. However, 
such state will be stabilized only in the low temperature 
and strong magnetic field region and will not change the 
essential feature studied in this paper. Therefore, we 
take the possible wave vectors of the order parameters 
as Q — Qo,2Qo,3Qq and 4Qo (equivalent to 0), where 
Qo = 2fc F . 

We evaluate the Helmholtz free energy by using the 
eigenvalue 6j and the unitary matrix C/kcj obtained by 
diagonalizing the mean field Hamiltonian K, +U M + V M + 
Vf 1 , where index j includes the degree of the spin free- 
dom. We determine the chemical potential p(p aa (Q,T)) 
from N — 1/4, where 



N= J2 <C\k x ,a)C{k x ,cx) > 

k x ,(T 

2N S 

= —Y 



cxp 1 £j - M 1 + 1 



T 



(11) 



The self-consistency condition (eq. (10)) is given by 
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where f(ej) is the Fermi function. 
(7) We obtain the free energy per site 



F( Paa (Q, T)) = 2pn - J2 lo S \ ex P 



2N S 



(8) 
V) 

(9) 



iVs 



T 



i^ m (Q,T)^(Q,T), 



U 
V 



m ]T e-^ a p aa (Q,T)p* a , a ,(Q,T) 



IP 

Yi 
u 2 



Q,<T,u' 



E e- 2lQa /w(Q,7>;, CT ,(Q,T), (13) 
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At T = it reduces to the ground state energy per site 



The electron density at the j site, n(j, T), its deviation 
from the mean value (1/2), 5(j, T), and the spin moment 
at j site, S z (j,T), are given by 



n(j,T) = ±Y,p^(Q,T)e iC » a = ± + 6(j,T), (15) 

Q,<7 



and 



S z (j,T) = ±-J2(Pn(Q,T)-p U (Q,T))e^ a . (16) 
Q 



When T = 0, n(j) = n(j,0), 5,(j) = &(j,0) and 
6(j)=6(j,Q). 




o. 2 - (t.t44> : 

o I — a* * x x x * H x x x 

1 2 V / t 3 4 5 

FIG. 1. 2S*(1) and 2S Z (2) as a function of V at H = 



0.5 




FIG. 2. 5 as a function of V at H = 



III. RESULTS AND DISCUSSIONS 
A. 2k F -SDW and 4k F -CDW 

In this subsection, we take Vi — and T = 0. It is 
known that the coexistent state of 2fcp-SDW and Akp- 
CDW are induced by U and V. |9). Figs. 1 and 2 are 
S z and (J as a function of V/t at £7/i = 5.0 and H = 0. 
For < U < 0.39, the antiferromagnetic order ((T,T, J.,.1), 
i.e., S z (l) = S z (2) = -S z (3) = -S z (4)) is stabilized and 
there is no charge order (8(1) = 8(2) = 6(3) = 6(4) = 0). 
The spin order of (T,T,j,J.) nas the wave vector of 2kp. 
For V/t > 0.39, the spin order becomes (T,0,|,0) (S z (l) = 
-S z (3), S z (2) = S z (4) = 0) and the charge order 
(6,-6,6,-6) coexists, where 5(1) = 5(3) = 5 and 5(2) = 
5(4) = -5. These orders, (T,0,|,0) and (5,-5,5,-5), 
mean 2fcp-SDW and 4/cp-CDW, respectively. These re- 
sults were obtained by Seo and Fukuyama. Q 

We show S z and 5 for U/t = 5.0 and V/t = ~ 5.0 
as a function of perpendicular field (h x = p^H/t) in 
Figs. 3 and 4. The antiferromagnetic state is grad- 
ually suppressed up to the critical field (h x =1.3, 1.3, 
1.5, 1.8 and 2.4 for V/t=Q, 1.0, 2.0, 2.5 and 5.0, respec- 
tively), as shown in Fig. 3, where S z (l) — —S z (3) and 
S z (2) — S z (4) — 0. Upon increasing h x the charge or- 
der decreases and becomes zero when V/t = 1.0, and 
the charge order becomes constant for h x > h x when 
V/t > 1.5, as shown in Fig. 4, where 6(1) — 5(3) = 6 and 
5(2) = 5(4) = —5. The /^-dependence of the amplitude 
of S z is almost the same as eq. (1) when we set H x as 
H x — th x /fiB, which is shown by solid lines in Fig. 3. 
From Figs. 3 and 4, we make the V-h x phase diagram, 
as shown in Fig. 5, where (0,0,0,0) means that the spin 
and charge orders do not exist. The dotted lines show 
the second order transition lines and the solid line show 



0.4 

N s /2 



e(h, p ^(q,o))= i J2 ^-TfEmW-ovti^o), 



N s ^ J U 

3 = 1 Q 0.3 

^ £ e-^/w(Q,0K, CT ,(Q,0) 6 
Q,o-,o-> o.2 

^ e'^ a p aa (Q,0)p:, a ,(Q,0). (14) 
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the first order transition lines. In the region of large h x 
and V, the state of 4/c F -CDW, (6,— 5,6,— 6), is stabilized. 



The h x -dependence of the spin order can be under- 
stood by the mean field solutions of the Ising model men- 
tioned in the introduction. |23|,[24| 



U/t = 5.0 




FIG. 3. 2S Z (1) as a function of h x . The solid line is given 
by eq. (1). 
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FIG. 4. 5 as a function of h x . 



3 V/t 4 



FIG. 5. V — h x Phase diagram. 



B. 2k F -SDW and 2k F -CDW 

Here, V% is introduced, which favors the coexistent 
state of 2fc F -SDW and 2fc F -CDW. (ll]Jl| We search the 
most stable solution at U/t = 5.0, V/t = 2.0, T = and 
H = by changing V% ■ These ^-dependences of the spin 
density and charge density are shown in Figs. 6 and 7. 
For < V 2 /t < 1.4, the coexistent state of 2fc F -SDW and 
4&F-CDW is stable, while the state changes to the co- 
existent state of 2fc F -SDW and 2fc F -CDW (( T ,T,l,i), i-e., 
S z (l) = -S z (4) < S z (2) = -5,(3)) and ((-6,5,6,-6), 
6(1) = 5(4) = -5 and 5(2) = 5(3) = 5) for V 2 /t > 1.4. 
This transition is a first order transition. 

The coexistent state of 2fc F -SDW and 2A: F -CDW ob- 
tained by Kobayashi et al. [0 is (TjTjIjI) anal ( — 5, 5, 
5, —5), which is different from the coexistent state of 
2/cp-SDW and 2fc F -CDW (( T ,T,l,i) and (-5, 6, 5, -6)). 
Since they limited the freedom of the order parameters, 
they could not find this coexistent state of 2/c F -SDW and 
2/cp-CDW (( T ,T,l,l) and (S, 5, 6, -6)). They also in- 
dicated that the coexistent state of 2/c F -SDW and 2fc F - 
CDW ((T,t,l,l) and (—5, 6, 5, —5)) is suppressed when 
the dimerization is not included. In fact, we could not 
find the coexistent state of 2/c F -SDW and 2/c F -CDW of 
(TjTiJ-)!) an< i (~5, 5, 5, —5) in the parameter region, 
(U/t = 5.0, V/t = 2.0 and V 2 /t < 2.0). However, as 
we show here, the coexistent state of 2/c F -SDW and 2fc F - 
CDW ((t,T,l,i) and (—5, 5, 5, —5)) is stabilized even in 
the absence of the dimerization. 

After we finished this study, we know the very recent 
study at T = @ and T ^ (29) in the absence of 
the magnetic field by Tomio and Suzumura. They cal- 
culate the mean field solution using the extended Hub- 
bard model with V and V 2 - They show that the state of 
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2fc F -SDW and 2fc F -CDW (( T ,T,l,i) and (-6, 6, 6, -6} 
coexists even if the dimerization is not included, too. 128 
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FIG. 6. 2S«(j) (j = 1 ■ ■ ■ 4) as a function of V 2 /t at H = 0. 
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FIG. 7. 5(j) (J = 1 • • -4) as a function of V 2 /t at # = 0. 



We analyze the coexistent state of 2fcp-SDW and 2kp- 
CDW under magnetic fields at < V 2 /t < 2.0. In 
the case of the perpendicular field (h x ), the state of 
(tjTjJ-u) and (-^,6,6,-6) becomes the one of (T,t,j,l) 
and (0,0,0,0) (6(1) = 6(2) = 6(3) = 6(4) = 0) at 0.75 < 
h x < 1.25, which is shown in Figs. 8 and 9. We can 
see that the coexistent state of 2/ep-SDW and 2fcp-CDW 
changes to the state of 2/cp-SDW. This transitions is the 
second order phase transition, since |5 Z (1)| and |5 Z (3)| 
(15*2(2)1 and 15*2(4)1) increase (decrease) gradually, and 
( T ,t,4,x) changes to (T,t,l,l) smoothly. When h x > 1.25, 
the 2/cp-SDW state disappears. For < V 2 /t < 1.4, 



as h x increases, the spin density in the coexistent state 
of 2fcp-SDW and 4fcp-CDW decrease and becomes zero 
at h% = 1.3,1.3,1.3 and 1.5 for V 2 /t = 1.4,1.0,0.5 and 
0, as shown in Fig. 10, where 5 2 (1) = — 5 2 (3) and 
5 2 (2) = 5 Z (4) = 0. When we set H° as H° = th c x /^ B , 
the /independence of S z is in agreement with eq. (1), 
which is shown by solid line in Fig. 10. The order 
of 4/cf-CDW becomes nonzero at the higher field for 
V 2 /t < 0.5, as shown in Fig. 11, where 6(1) = 5(3) = 5 
and 6(2) = 5(4) = —5. We show the V 2 -h x phase in Fig. 
12, where the solid line is for the first order transition 
and the dotted lines are for the second order transitions. 




FIG. 8. 2S z (j) (j = 1 • • -4) as a function of h x . 
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FIG. 9. 5(1) and 6(2) as a function of h x . 
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FIG. 12. V — h x phase diagram. 



C. Finite Temperature 




FIG. 11. 8 as a function of h x 



We show S z and 5 as a function of T at 7? = 0. At 
[//t = 5.0 and Vz/t = 0, 5, and 5 as a function of T for 
various values of V are shown in Figs 13 and 14, where 
5,(1, T) = S,(2,T) = -5,(3, T) = -S Z (A,T) for V/t = 
and S g (l,T) = —S z (3, T),S Z (2, T) = 5,(4, T) = 0, 
6(1, T) = 5(3, T) = 5(T) and 5(2, T) = 5(4, T) = -5(T) 
for V/t — 0.5 ~ 5.0. As V increases, the critical temper- 
atures (Tsdw and Tcdw) a t which S z and 6 become zero 
increase. The T-dependences of n for V/t — 3.0, 4.0 and 
5.0 have a dip at Tsdw for V/t = 3.0, 4.0 and 5.0, which 
can be seen in Fig. 14. We show the V-T phase diagram 
in Fig. 15, where the solid and dotted lines are for the 
first and the second order transitions, respectively. Note 
that Tcdw > ?sdw for V/t > 2.0, that is, the charge or- 
der remains even if the spin order disappears under high 
temperatures. 

As we have shown above, the coexistent state of 2kp- 
SDW and 2fc F -CDW (( T ,T,4u) and (-5,5,5,-5)) are sta- 
bilized for V 2 /t > 1.4 at T/t = (see Figs. 6 and 7). 
At T/t ^ 0, for example, for Vi/t = 1.6 as T increases, 
the magnitudes of S z and n decreases and the second or- 
der transition to the state of 2fcp-SDW ((1\1\Jo!)) occurs, 
as shown in Figs. 16 and 17. Namely, T CD w < ?sdw 
in the coexistent state of 2fcp-SDW and 2fcp-CDW. For 
Vi/t < 1.4, Tcdw m the coexistent state of 2fcp-SDW and 
4fcp-CDW is the same as Tsdw, as shown in Figs. 18 and 
19, where S Z (1,T) = -5,(3, T), 5,(2, T) = 5,(4,T) = 0, 
5(1, T) = 5(3, T) = 5(T) and 5(2, T) = 5(4, T) = -5(T). 
We show the V2-T phase diagram in Fig. 20, where the 
solid and dotted lines are for the first and the second 
order transition lines, respectively. 

Our phase diagrams (V-T and V2-T) are in good agree- 
ment with those obtained by Tomio and Suzumura. |29[] 



6 



\\ U/t=5.0 






■ v/t=o 


- -□- 


■ V/t = 1 .0 


- -o- 


■ V/t=2.0 


- -X- 


■ V/t=3.0 


- - +- 


- V/U4.0 


- - A- 


■ V/t=5.0 



T/t 

FIG. 13. 2S Z (1,T) as a function of T/t. 
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FIG. 15. V — T Phase diagram. 
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FIG. 14. S(T) as a function of T/t. 



FIG. 16. 2S Z (J, T) (j = 1 • • • 4) as a function of T/t. 
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FIG. 17. 5(1, T) and 5(2, T) as a function of T/t. 

FIG. 19. 8{T) as a function of T/t. 
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FIG. 18. 2S Z (1,T) as a function of T/t. 



FIG. 20. V 2 - T Phase diagram. 



IV. CONCLUSIONS 

We have shown the different dependences of the 
magnetic field and temperatures on the 2fc F -SDW-2A:F- 
CDW coexistent state and 2fc F -SDW-4fc F -CDW coexist- 
ing state. 

When V2 is absent or small, the coexistent state of 
2fc F -SDW and 4fc F -CDW is stable. Then the charge or- 
der survives at high fields in the case of the perpendic- 
ular magnetic field to the easy axis when V/t > 2.0 . 
This situation may be realized in (TMTTF) 2 X, where 
4fc F -CDW is expected to be observed in X-ray scattering 
measurements under high fields. 
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On the other hand, if V2 is large enough, the coexistent 
state of 2fc F -SDW and 2fc F -CDW is stabilized. When the 
magnetic field is applied perpendicular along the easy 
axis, the coexistent state of 2/cp-SDW and 2fcp-CDW 
changes to 2/cp-SDW state. This transition is the second 
order phase transition, and (|,T,l,i) becomes (T,f,j,.l)- 
This transition may be observed in the angle dependence 
of satelite peak positions of NMR in (TMTSF) 2 PF 6 . 

The critial temperature of the charge order (Tcdw) is 
higher than that of the spin order (Tsdw) m the coex- 
istent phase of 2A: F -SDW and 4fc F -CDW when V/t > 
2.0. This is in good agreement with Tsdw ~ 13 K 
and Tcdw ~ 100 K observed in X-ray scattering in 
(TMTTF) 2 Br. Since T CD w is convergent to T SDW for 
V/t < 2.0, we expect that the critical temperatures of 
2A: F -SDW and 4/cp-CDW will become the same when the 
pressure is applied to (TMTTF) 2 Br. 

In the coexistent phase of 2/cf-SDW and 2fcF-CDW, 
we find that Tcdw < ^SDWj which is consistent with the 
X-ray scattering measurement in (TMTSF^PFg, where 
^cdw — ^sdw- M 
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